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Energetics of the fee (7) - bee (a) lattice transformation by the Bain tetragonal deformation is 
calculated for both magnetically ordered and paramagnetic (disordered local moment) states of iron. 
The first-principle computational results manifest a relevance of the magnetic order in a scenario of 
the 7 - a transition and reveal a special role of the Curie temperature of a-Fe, Tc, where a character 
of the transformation is changed. At a cooling down to the temperatures T < Tc one can expect 
that the transformation is developed as a lattice instability whereas for T > Tc it follows a standard 
mechanism of creation and growth of an embryo of the new phase. It explains a closeness of Tc to 
the temperature of start of the martensitic transformation, Ms- 

PACS numbers: 61.50.Ks, 64.70.kd, 75.30.Et, 75.50.Bb, 71.20.Be 



Deeper understanding of mechanisms of polymorphous 
J - a transformation in iron and its alloys is of fundamen- 
tal importance for both metallurgical technologies [H, Q 
and for a general theory of phase transitions in solids 
[1) 0, 01 • Despite numerous investigations an issue of a 
mechanism of a new phase nucleation in the course of 7 
- a transformation remains open (see, e.g., discussion in 
Ref. i). 

It is well known fl*! that a character and rate of the 
transformation in iron and iron-based alloys is changed 
drastically below some temperature Ms, namely, at T < 
Ms it occurs by a fast cooperative shear deformation of 
atoms (martensitic mechanism) whereas for higher tem- 
peratures (but lower than the temperature of phase equi- 
librium Tj-a) the transformation develop much slower 
and a formation and growth of grains of a-phase is ob- 
served. In pure iron the temperature of start of the 
martensitic transformation Ms ~ 1020iir 01 which is 
only 30 K lower than the Curie temperature^ of a-Fe, Tc, 
in accordance with an old idea of Zener |S] suggesting 
that Ms ^ Tc- It is know also that in the system Fe-C a 
process of the j - a transformation becomes much faster 
at a cooling down below Tc jSl ■ It is commonly accepted 
now that magnetic degrees of freedom play a crucial role 
in the phase equilibrium in iron 

[IS El, Hi [II- However, 
a mechanism of their effect on kinetics of the transforma- 
tion remains unclear. 

Several ways to transform the crystal lattice from fee 
(7) to bcc (a) structures have been suggested, among 



them schemes of Bain [14] (a tetragonal deformation 
along the < 001 > axis) and of Kurdumov-Zaks [1] (two 
shear deformations) are the best known. To determine 
the deformation path and estimate energy barriers calcu- 
lations of total energy in a configurational space of lattice 
deformations are required. For the Kurdumov-Zaks path 
this is a rather cumbersome and computationally expen- 
sive problem. Therefore first-principle calculations of en- 
ergetics of polymorphous transformation in iron were car- 
ried out, up to now, only for the Bain transformation 
path [l3,[lB,[ll,[l3,lli,[il[20|- It turned out that ener- 
getics of the transformation is essentially dependent on 
magnetic structure of iron. In particular, ferromagnetic 
(FM) state of the fee iron is unstable with respect to 
tetragonal deformation, with energy minima at c/a ratio 
equal 1 (bcc) or c/a ~ 1.5 (fct structure) 0, [3, 
Antifcrromagnetic (AFM) state of the fee iron has lower 
energy than the FM one and show a monotonous increase 
of energy along 7 — > a transformation path 15, 16| . 



Since the magnetic ground state is different for a- and 
7-Fe, an investigation of energetics of the transformation 
path for a given coUinear magnetic structure [13, O, 

[i3,[Ii,[3 

is not enough to describe properly the a — 7 
transition. As was shown in Refs. [2ll, |22, l23| the true 
magnetic ground state of fee Fe is rather complicated 
and, in general, noncoUinear (depending on the lattice 
constant). Evolution of magnetic state of Fe along the 
Bain path was studied in Ref. 20 by the TB-LMTO 
method in atomic sphere approximation (ASA) [13]. It 



2 



was demonstrated that, for an essential part of the Bain 
path, noncoUinear magnetic structures take place which 
are replaced by FM ordering for c/a smaller than some 
critical value. 

iSlii, 



The calculations 15 



19| have been carried out 
for a fixed value of volume per atom or lattice constant 
a that corresponds to so-called epitaxial Bain path de- 
scribing the transformation for iron films on a substrate. 
The latter leads to some restrictions on geometry of the 
transformation which takes place at low temperatures in 
magnetically ordered state. The results of Ref. allow 
us to predict a type of magnetic structures which can be 
realized in the epitaxial iron films. At the same time, the 
"f — a transformation in the bulk occurs at high tempera- 
tures when the magnetic state is disordered. To describe 
this situation, calculations for paramagnetic iron are re- 
quired. 

In this work we analyse in detail the Bain deforma- 
tion path for both noncoUinear and paramagnetic states 
using the methods of spin spirals (SS) and disordered lo- 
cal moments (DLM), respectively. As a result, we clarify 
the reasons for essential differences in mechanisms of the 
polymorphous transformation below and above the Curie 
temperature Tc- 

The calculations have been carried out using the VASP 
(Vienna Ab-initio Simulation Package) [H, [2^, with 
first-principle pseudopotentials constructed by the pro- 
jected augmented waves (PAW) [2^, and the generalized 
gradient approximation (GGA) for the density functional 
in a form (29l | with the parametrization (30l | (see also com- 
ments in Ref. ^ concerning use of the GGA for non- 
coUinear case). The PAW potential with energy mesh 
cutoff 530 eV, and uniform fc^oint 12x12x12 mesh in 
the Monkhorst-Park scheme [S^ were used. 

To describe the noncoUinear magnetic state the model 
of flat SS 3l| has been used with the magnetization rota- 
tion around wave vector of SS, q which chosen along the 
axis of the Bain tetragonal deformation < 001 >. The 
values q ~ and (7 = 0.5 (in units of 27r/c) correspond to 
the FM and AFM states, respectively. To calculate the 
energy of Fe with SS magnetic structure we employed the 
PAW formalism described in Ref. |34| and implemented 



in the VASP code. We optimized the volume per atom 
for given values of q and tetragonal deformation. 

Paramagnetic state of iron was modeled by the disor- 
dered local moment (DLM) method [sH]. To this aim, we 
used 27-atom supercell (Fig. 1) with a given random dis- 
tribution of magnetic moments (with zero total magnetic 
moment) which was kept fixed at the self-consistency pro- 
cess. The latter was provided by using constrained den- 
sity functional approach As well as for the SS case 
we optimized the volume per atom for a given tetragonal 
deformation. 

Computational results for SS are shown in Fig. 2. One 
can see that the minimum of the total energy is reached 
for bcc FM (curve 2) and fct AFM (curve 8), with a tran- 




FIG. 1: Supecell used in the DLM calculation. Arrows show 
the magnetic moments. 
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FIG. 2: Dependence of the total energy per atom for SS 
on the wave vector q ||< 001 > for different values of 
tetragonal deformation. Curves 1 to 9 correspond to c/a = 
0.9; 1.0; 1.2; 1.3; 1.35; V2; 1.5; 1.6; 1.7, respectively; c/a = 1 for 
bcc and c/a — \/2 for fee structures. 



sition between these two states via noncoUinear magnetic 
structure with the wave vectors q ~ 0.25 — 0.35. The en- 
ergy difference between the FM and AFM states is small 
for 1.3 < c/a < 1.5 and increases, almost in order of 
magnitude, for c/a < 1.3. This agrees qualitatively with 
the previous resuUs [20] of the TB-LMTO-ASA calcula- 
tions that the tetragonal deformation induces a magnetic 
transition. 

Similar to Refs. 2^ 31, 3^ we found that the energy 
minimum for the fee Fe corresponds to g ~ 0.3 (curve 
6 in Fig. 2). Optimization of the volume makes the 
q-dependence of the total energy rather flat, thus, the 
energy at g = 0.3 is only 6 meV/atom lower than for the 
AFM fee state. The FM fee structure {q = 0) has the 
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energy in 54 meV/atom higher than for q = 0.3. It is 
characterized by much larger values of magnetic moment 
and volume per atom (11.86 A'^ and 2.3 /xb at g = vs 
10.74 and 1.4 /z^ at g = 0.3). One can see in Fig. 
2 that the FM state of fee iron is unstable with respect 
to the tetragonal deformation; if one keeps the FM order 
the fee iron will reconstruct spontaneously to bee or fct 
phase, without essential changes of volume and magnetic 
moment per atom. 

Fig. 3 presents the total energy along the Bain path for 
optimized (with a given c/a ratio) SS magnetic structure, 
together with the results for the coUinear (FM and AFM) 
case, for the paramagnetic case (DLM) and for partially 
disordered magnetic structure (DLM0.5). The results for 
the FM and AFM magnetic structures are close to those 
obtained earlier by full-potential FLAPW method [3]. 
Taking into account noncoUinear magnetic configurations 
allows us to describe a continuous transition between FM 
and AFM states and to estimate the energy barrier Ei, 
resulting from the Bain deformation. We have found 
Eb ~ 20 meV/atom for fcc-bcc and Ej, ~ 40 meV/atom 
for fct-bcc transitions. The corresponding values found 
in Ref. [l^ are 48 meV/atom and 63 meV/atom, respec- 
tively which is essentially larger than our data. The dif- 
ference is mainly due to optimization of volume per atom 
used in our calculations (in Ref. fixed values of the 
volume were used). 
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FIG. 3: Variation of the total energy per atom along the Bain 
deformation path for different magnetic states. FM (empty 
triangles) and AFM (solid triangles) label colhnear ferromag- 
netic and antiferromagnetic structures, SS (solid circles) - 
spin-spiral state, DLM (empty circles) - disordered local mo- 
ments, DLM0.5 (empty diamonds) - DLM state with the total 
magnetic moment equal to half of that for the FM state. 

Comparing Figs. 2 and 3 one can see that the barrier 
position corresponds to the value of tetragonal deforma- 
tion where magnetic state switches from SS {q ~ 0.3) to 
FM {q = 0). Interestingly, the critical value of the Bain 



deformation found for the SS state turns out to be close 
to the crossing point of the corresponding curves for FM 
and AFM configurations. 

For the paramagnetic (DLM) state the energy has a 
minimum corresponding to the fee structure and is close 
to the energy of SS and AFM states for 1.3 < c/a < ^2 
(Fig. 3). For c/a < 1.3 the energy grows monotonously 
with the c/a decrease reaching the maximum for the bcc 
state (c/a = 1). This maximum is much higher than the 
energy for FM bcc and DLM fee configurations (by 220 
and 110 meV/atom, respectively). Thus, reconstruction 
of crystal lattice of Fe from fee to bcc by the Bain de- 
formation without change of magnetic configuration does 
not lead to any energy gain. To model the state of iron 
at finite temperatures below Tc we have performed also 
calculations for the DLM state with total magnetization 
M — O.bMmax where Mmax is the magnetization in FM 
state for given c/a ratio. The results are not essentially 
different from the paramagnetic case except for region 
close to bcc structure (compare the curves DLM and 
DLM0.5 in Fig. 3). So, the energetics of the Bain trans- 
formation path changes drastically when the magnetic 
state becomes close to the ferromagnetic. 

It is not surprising, of course, that the energy of bcc 
Fe is much higher in the DLM state than in ferromag- 
netic one, in light of a common opinion (originated from 
a seminal work by Zener (37| ) that it is ferromagnetism 
that stabilizes a-Fe. According to Ref. [H a strong 
enough short-range magnetic order survives in iron above 
Tc ~ 1043 K which is probably essential to explain a 
stability of bcc phase in some temperature region above 
Tc [Hi [s^l ■ To consider this problem quantitatively one 
has to be able to calculate phonon and magnetic con- 
tributions to entropy of different phases which is rather 
complicated. However, since Tc is much higher than typ- 
ical phonon energies the lattice vibration entropy in both 
phases are approximately constant in the temperature re- 
gion under consideration and therefore the phonon con- 
tribution to the difference of free energies between a- and 
7-phases is a smooth function linearly dependent on the 
temperature. At the same time, magnetic contribution to 
the total energy is strongly dependent on the magnetic 
configuration and thus on the temperature, as follows 
from the presented results. 

When cooling down iron in a temperature interval 
Tc < T < T-y_Q, the bcc structure becomes preferable 
due to entropy contribution in free energy. Since the dif- 
ference between free energies of the bcc and fee phases 
is zero at T-y-a and changes slowly with the tempera- 
ture, the moving force of the phase transition is relatively 
weak. As a result, the 7 — 0; transformation develops 
in this situation by the classical nucleation and growth 
mechanism. 

When decreasing the temperature below Tc the FM 
state of bcc Fe arises which has the energy much lower 
than that of paramagnetic fee state (by 110 meV/atom). 
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At the same time, lattice reconstruction from paramag- 
netic fee to FM bee state requires to overcome a rather 
low barrier which height (~ 20 meV/atom, ot ^ 250 
K/atom) is small in comparison with the temperature 
T ^ Ms- Thus one can expect that cooling down to 
T < Tc will initiate martensitic mechanism of 7 — a 
transformation, via a development of lattice instability 
I39II . This conclusion differs essentially from that of Ref . 
I12I where use the less accurate LMTO method with model 
Stoner parameters has resulted in three times larger bar- 
rier height than in our calculations. 

Thus, Tc of a-Fe plays a role of a special point where 
kinetics of the transformation is changed, due to a cru- 
cial role of magnetic degrees of freedom in energetics of 
iron which is a main conclusion from our computational 
results. 
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